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OBSERVATIONAL CONSTRAINTS 
ON COSMOLOGICAL MODELS WITH CHAPLYGIN GAS 
AND QUADRATIC EQUATION OF STATE 
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^ Tver state university 
170002, Sadovyj per. 35, Tver, iJwssicQ 

Observational manifestations of accelerated expansion of the universe, in particular, recent 
data for Type la supernovae, baryon acoustic oscillations, for the Hubble parameter H{z) 
and cosmic microwave background constraints are described with different cosmological mod¬ 
els. We compare the ACDM, the models with generalized and modified Chaplygin gas and 
the model with quadratic equation of state. For these models we estimate optimal model 
parameters and their permissible errors with different approaches to calculation of sound 
horizon scale rs{zd). Among the considered models the best value of is achieved for the 
model with quadratic equation of state, but it has 2 additional parameters in comparison 
with the ACDM and therefore is not favored by the Akaike information criterion. 


I. INTRODUCTION 


Observations 0,0 of Type la supernovae demonstrated accelerated expansion of our universe. 
Further investigations of supernovae i,a, baryon acoustic oscillations 


background measurements |7h1C1||. estimations [lll-l25l| of the Hubble parameter H{z) for different 


cosmic microwave 


redshifts z conhrmed accelerated growth of the cosmological scale factor a(t) at late stage of its 
evolution. 

For Type la supernovae we can measure their redshifts 2 and luminosity distances Dp, so these 
objects may be used as standard candles 0-0. 

Baryon acoustic oscillations (BAO) are observed as a peak in the correlation function of the 
galaxy distribution at the comoving sound horizon scale rs{zd) 0,0, corresponding to the end of 
the drag era, when baryons became decoupled and acoustic waves propagation was ended. This 
effect has various observational manifestations E-[l(il . Mi, in particular, one can estimate the 
Hubble parameter H{z) for dehnite redshifts |16l425j| (details are in Sect. HI)) . 

The mentioned recent observations of Type la supernovae, BAO and H{z) essentially restrict 
possible cosmological theories and models. To satisfy these observations all models are to describe 
accelerated expansion of the universe with dehnite parameters @-[i3,&[3i. 

The standard Einstein gravity with A = 0 predicts deceleration of the expanding universe: 
a"{t) < 0. So to explain observed accelerated expansion, we are to modify this theory. The most 
simple (and most popular) modihcation is the ACDM, including dark energy corresponding to 
A 7 ^ 0 and cold dark matter in addition to dehcient visible matter. This model with appropriate 
parameters 0@-[l0 successfully describes practically all observational data, in Sect. IHII we apply 
this model to describe the updated recent observations of Type la supernovae, BAO effects and 
H{z) estimates. In this paper we use the notation ACDM for the model with an arbitrary spatial 
curvature Dfc- 

One should note that there are some problems in the ACDM model, in particular, ambiguous 
nature of dark matter and dark energy, the problem of hne tuning for the observed value of A and 
the coincidence problem for surprising proximity Da and nowadays 35|, . 
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Therefore cosmologists suggested a lot of alternative models with different equations of state, 
scalar fields, with f{R) Lagrangian, additional space dimensions and many others In this 

paper we consider in detail two models with nontrivial equations of state describing both dark 
matter and dark energy: the model with modified Chaplygin gas (MCG) [40-47] in Sect. m and 
the model with quadratic equation of state 4^5^ in Sect. |Vl 


II. OBSERVATIONAL DATA 


In this paper we use the Union2.1 compilation of Type la supernovae (SNe la) observational 
data. This table includes redshifts z = Zi and distance moduli //j = with errors cjj for 
^SN = 580 supernovae. The distance modulus //j = h[Dl) = Slogig (.D^^/lOpc) is logarithm of 
the luminosity distance [ij, IH, 35 j: 


Dl{z) = 


c{l + z) 
Ha 


Sk[H, 


dz 

W) 


( 1 ) 


Here 

{ sinh(xv^)/V^, > 0, 

X, Qk = 0, 

sin (xA/|Ofc|)/\/|Hfc|, Ofc < 0; 


redshift z and the Hubble parameter H{z) are connected with the scale factor a{t): 


a{t) 


ao 

1 + z’ 


H(z) 


a(i) . 

a(t)’ 


( 2 ) 


k is the sign of curvature, present time fraction, og = a(to) and Hq = H{to) 

are the current values of a and H. 

For any cosmological model we fix its model parameters pi,p 2 , ■ ■ ■, calculate dependence a{t), 
the integral ([T]) and this model predicts theoretical values for luminosity distance ([T]) (for given 
z), or for modulus. To compare these theoretical values with the observational data Zi and 
from the table we use the function 


Nsn 

x|iv(Pi,P2,...) = min XI A/i* =F7g,pi,...) -/if *. (3) 

ij=i 

Here Csn is the SN-by-SN covariance matrix [^, representing systematic errors. 

In the sum m marginalization over the Hubble constant Hq is assumed, because we have to 
take into account model de^ndence of the moduli /if®. Unlike observed apparent magnitudes 
mf® the values /if® in Ref. are estimated as 

^obs ^ _ M + axi- (3c + 6P. (4) 


This formula includes the SN la absolute magnitude M and corrections connected with deviations 
from mean values of lightcurve shape (xi), SN la color (c) and mass of a host galaxy (the factor 
P). The parameters M, a, (5 and <5 are considered in Ref. as nuisance parameters, they are 
fitted simultaneously with the cosmological parameters in the flat ACDM model. 

Thus we have a model dependent additive term in Eq. (jl]) for the Union2.1 values /i°^® with 
concealed dependence on the Hubble constant Hq and other model parameters. In particular, one 
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can find only that the best fit value for the absolute magnitude M = —19.321 ± 0.03 is obtained 
in Ref. for = 0.7, where h = 77o/100 kms~^Mpc~^. 

To suppress this dependence many authors 54457t| suppose that values from any SN la 
survey have a systematic error depending on Hq and marginalize the sum ([T]) over the Hubble 
constant Hq. They use the fact, that for the most popular models theoretical value of the luminosity 
distance ({!]) depends on Hq as ~ H^ , so the distance modulus has the additive term 
—Slog^o-^o- In Ref. [54] this term is separated as ^q in the form 


h{Dl) = 5 logio 1*0 = 42.384 - 5 log^o h] 

if we denote d* = A^Uj — /Uq, d = {di,... ,dsN), the minimum of the sum ([3]) over Hq (or over /io) 
will take the form 


XsNiPu---) 




id=i 


id=i 


(5) 


In this paper for all models we use the marginalized function ([5]) Xsn 1° describe the supernovae 
la data [^. 

This approach with separation of the Hubble constant Hq among other model parameters can 
not be applied to H{z) and BAO observational data, because observed values have different de¬ 
pendence on Hq. 

To describe the BAO data we calculate the distance 0,i-0 


Dv{z) 


czDli 


-,1/3 


{l + zyH{z) 


( 6 ) 


and two measured values 


dz{z) = 


rsizd) 

Dv{z)^ 


A/ \ TdoVIIm p, / N 
A{z) = - Dv{z), 


cz 


(7) 


which are usually considered as observational manifestations of baryon acoustic oscillations 
Here = ^TTGp{tQ)/H q is the present time fraction of matter with density p. The value rs{zd) 
in Eq. (0) is sound horizon size at the end of the drag era Zd- 



( 8 ) 


To estimate this important parameter different anthors 
proaches and snggested different fitting formulas for r*. 
of rs(zd) = Td and r^h are shown: 


17h 33I| used theoretical or statistical ap- 


In table H] the following recent estimations 


Refs 


\1,M 

[10] 


[33] 

[22] 

[29] 

m 

[2^ 

\27,28,M 

m 

Tsizd) 

147.4 

147.49 

147.6 

148.6 

148.69 

149.28 

152.40 

152.76 

153.19 

153.2 

153.5 

Td ■ h 

98.79 

99.26 

100.09 

99.78 

99.62 

104.49 

105.15 

106.93 

107.23 

103.6 

104.69 


TABLE I: Recent estimations of rs{zd) and rs{zd) ■ h (Mpc). 


In the most of cited papers in table U the values rs{zd) were considered as fiducial ones for 
calculating Dv{z), H{z) and other parameters. So these results sufficiently depend on rs{zd), in 
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particular, estimations of H{z) from the BAO data 1^25| are defined with the factor the 

corresponding factor rd/r^d^ takes place for calculated values Dy{z) or Da{z) = Dl{z)/{ 1 + z)^. 

In this paper we use two different approaches to calculate the sound horizon scale rs{zd)- But 
previously one should mention the simplest method, applied in Ref. (s^ . where the arithmetic 
average of the values in table H] (with their multiplicity) 


r^izd) = 150.69 ±2.45 Mpc 


(9) 


was used as the basic value. This value is independent on Hq, hence the observational parameter 
dz{z) in Eq. ([7|) appears to be Hubble dependent (though the formula ([8|) predicts ~ and 
Hubble free dz)- One may conclude, that h dependence of dz is the drawback of this approach, 
so in this paper we consider the fixed value ([9]) only in section IHII to emphasize advantages of 
other methods. 

More appropriate procedures to calculate rs{zd) include different fitting formulae B 0,0 [ 53 . 
In this paper we use the numerically calibrated approximation from Ref. [TJ 


= 


55.154 exp [72.3(H^/i2 + 0.0006)^ 


Mpc 


( 10 ) 


(H™/i2)0-25351(J]^/j2)0.12807 

as the basic formula. The resulting h dependence in Eq. (1101) (for a reasonable neutrino contribution 
with '>^ 1 ' ^ 0-23 eV 0 ) is Td ^ h 0.7632^ ji- -g close to the true variant rd ^ h The 

dependence on in Eq. (|10p is well fitted for ACDM-like models, however for the models with 
Chaplygin gas and with quadratic EoS, considered below, 11^ is not a basic model parameter. 
The value in these models should be estimated in a special way, so an additional uncertainty 
appears in this approach. 

Thus, an alternative simple fitting formula 


rsizd) = 


iXd ■ h) fid 
h ’ 


{rd ■ h)fid = 104.57 Mpc. 


( 11 ) 


with true h dependence may be suggested. Here the value (rd ■ h)fid = 104.57 ± 1.44 Mpc was 
chosen as the best fit for the ACDM model. This procedure is described in the next section and 
illustrated in figure [H 

The parameter (r^ • h)fid for the expression (|lip plays the same role as the baryonic fraction 
for the formula (1101) . in both cases we do not consider and as free model parameters for all 
models, but fix their optimal (fiducial) values after description of the simplest ACDM model. The 
best ACDM fit for in Eq. (fTOP (see figure [I|) is 


Qh = 0.044 ± 0.004. 


( 12 ) 


To take into account all available BAO data 0-0 for parameters ([7|), we consider in this 
paper Nbao = 17 data points for dz{z) (10 additional points in comparison with the table in our 


paper 


39l |) and 7 data points for A{z) presented in the table UTl 


Measurements of dz{z) and A{z) from Refs. 2^, l29j in table [H] are not independent. So the yX 
function for the values ([7]) is 

xIao{pi.P2,---) = (Ad)^C7'Ad± (AA)^C^iAA, 


The elements of covariance matrices 


and C A = 


-1 — lUAi 


cfa = 30124, 4 = -17227, = 86977, 

c'Xn = 24532.1, 4^4 = -25137.7, = 12099.1, 

cf 4 j 5 = -64783.9, c^i 5 = 128837.6; = 1040.3, 

= 336.8, 4^4 = 3720.3, = -1551.9, 


Ad = dzizi) - df. (13) 

in Eq. (fl^ are BS0: 

4^4 = 134598.4, 

4i4 = -807.5, 


-BK 


= 2914.9. 
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z 

dziz) 

CTd 

Aiz) 

(TA 

Refs 

Survey 

0.106 

0.336 

0.015 

0.526 

0.028 

[8,^ 

6 dFGS 

0.15 

0.2232 

0.0084 

- 

- 

[3^ 

SDSS DR7 

0.20 

0.1905 

0.0061 

0.488 

0.016 

[26,^ 

SDSS DR7 

0.275 

0.1390 

0.0037 

- 

- 

\M 

SDSS DR7 

0.278 

0.1394 

0.0049 

- 

- 

[27] 

SDSS DR7 

0.314 

0.1239 

0.0033 

- 

- 

[29] 

SDSS LRG 

0.32 

0.1181 

0.0026 

- 

- 

[22] 

BOSS DRll 

0.35 

0.1097 

0.0036 

0.484 

0.016 

[26,^ 

SDSS DR7 

0.35 

0.1126 

0.0022 

- 

- 

[30] 

SDSS DR7 

0.35 

0.1161 

0.0146 

- 

- 

\m 

SDSS DR7 

0.44 

0.0916 

0.0071 

0.474 

0.034 

[29] 

WiggleZ 

0.57 

0.0739 

0.0043 

0.436 

0.017 

[20] 

SDSS DR9 

0.57 

0.0726 

0.0014 

- 

- 

[22] 

SDSS DRll 

0.60 

0.0726 

0.0034 

0.442 

0.020 

[29] 

WiggleZ 

0.73 

0.0592 

0.0032 

0.424 

0.021 

[29] 

WiggleZ 

2.34 

0.0320 

0.0021 

- 

- 

[2a 

BOSS DRll 

2.36 

0.0329 

0.0017 

- 

- 

m 

BOSS DRll 


TABLE II; Values of dz{z) = Tsizd)/Dv{z) and A{z) ([T]) with errors and references 


Here Cij = Cji, the remaining matrix elements are ca = l/cr?, dj = 0, i ^ j. 

In the values Cd in table El we took into account correlation between estimations of dz{z) and 
H{z) (table Uni) for z = 0.35, 0.57, 2.34, 2.36 in Refs. [iMl|,0,[^|. 


Measurements of the Hubble parameter H(z) for different redshifts z with 38 data points 11142 f 
are presented in table III^ These values H{z) were calculated with two methods: 1) differential age 
approach in Refs. Ill4l5l| with evaluation of the age difference dt for galaxies with close redshifts 
dz and the formula 


Hiz) = 


1 da 


1 dz 


a{t) dt 1 + z dt' 

2 ) measurement U25| of the BAO peak in the correlation function in line-of-sight directions at 
a redshift separation Az = rs{zd) H{z)/c. 

3 "^ of H{z) essentially depend on a fiducial value 


For the latter method estimations 


and have the factor /vd, as was mentioned above. In particular, the result in Ref. 25] is 


H{z = 2.34) = (222 ± 7) 


km 147.4 Mpc 


Mpc TsiZd) 

In table uni this factor is taken into account only for the errors an from the papers 0, [3- 
0, 24, 25], where fiducial values essentially differ from the average Q. The estimations for 


H{z) in table Ell are the same as in the correspondent sources. 

To compare the H(z) data in table Ell with = 38 data points with model predictions we 
use the function 


Xh{pi,P2,---) = 




i=l 


a 




similar to the function ([3]) for the SN la observational data from Ref. [^. 


(14) 
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z 

H{z) 

cth 

Refs 

z 

Hiz) 

cth 

Refs 

0.070 

69 

19.6 


0.570 

96.8 

3.4 

[22] 

0.090 

69 

12 

[1^ 

0.593 

104 

13 


0.120 

68.6 

26.2 

\u\ 

0.600 

87.9 

6.1 

[17] 

0.170 

83 

8 

[1^ 

0.680 

92 

8 

\m 

0.179 

75 

4 

\m 

0.730 

97.3 

7.0 

[17] 

0.199 

75 

5 

[1^ 

0.781 

105 

12 

m 

0.200 

72.9 

29.6 

M 

0.875 

125 

17 

[1^ 

0.240 

79.69 

2.99 

[1^ 

0.880 

90 

40 

[12] 

0.270 

77 

14 

[1^ 

0.900 

117 

23 

[1^ 

0.280 

88.8 

36.6 

\M 

1.037 

154 

20 

\m 

0.300 

81.7 

6.22 

\M 

1.300 

168 

17 

[1^ 

0.340 

83.8 

3.66 

[1^ 

1.363 

160 

33.6 

[la 

0.350 

82.7 

9.1 

m 

1.430 

177 

18 

[1^ 

0.352 

83 

14 

[1^ 

1.530 

140 

14 

[1^ 

0.400 

95 

17 

[1^ 

1.750 

202 

40 

[1^ 

0.430 

86.45 

3.97 

\m 

1.965 

186.5 

50.4 

[la 

0.440 

82.6 

7.8 

[17] 

2.300 

224 

8.6 

m 

0.480 

97 

62 

[12] 

2.340 

222 

8.5 

[2^ 

0.570 

87.6 

7.8 

[20] 

2.360 

226 

9.3 

m 


TABLE III; Values of the Hubble parameter H{z) with errors an from Refs. 


11H25| 


We mentioned above that the observed values in tabies mm have different dependence on 
h = Hq/100 kms“^Mpc“^. In particular, if we use the fitting formula (fTO]) for rs{zd), the parameter 
d, © has rather weak h dependence {dz ~ /jO-2368^ because Dy ^ h ^); for the formula (HH) dz 
is Hubble free in accordance with Eq. ©. The value A(z) is also Hubble free, but it depends on 
Clm- The values H{z) in table IHlI are naturally proportional to Hq. So estimations in tables HIl IHII 
may be model dependent, but we can assume that different authors use different methods and 
produce possible systematic errors for dz{z) and H{z) in different directions. One can suppose 
mean systematic errors to be close to zero. 

On the other hand, if any form of marginalization over Hq for BAO and H(z) data is made [si- 


57l |. the obtained results will have an additional error, because a model can successfully describe 


all SN la, BAO and H(z) data, but with 3 essentially different intrinsic values of Hiy Under these 
arguments we make the marginalization procedure ([5]) over Hq only for SN la data [3|], but not for 
BAO and H{z) data from tables HH IHII 


III. ACDM MODEL 

In the ACDM and other models in this paper the Einstein equations 

= SttGTI; - (15) 

describe dynamics of the universe. Here G^ = Tj^ = diag{—p,p,p,p). 

In the ACDM model baryonic and dark matter may be described as one component of dust-like 
matter with density p = Pc = Pb + Pdrm so we suppose p = 0 in T^. In models with Chaplygin gas 
fSect. irvTl and with quadratic equation of state (Sect.|V]) we suppose that an additional component 
of matter describes both dark matter and dark energy and gives some contribution pg in the total 
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density: 


P — Pc + Pg + Pr- 


(16) 


The fraction of relativistic matter (radiation and neutrinos) is close to zero for observable values 
z < 2.36, so below we suppose pr = 0. 

For the Robertson-Walker metric with the curvature sign k 


ds^ = —dt^ + a^{t) (1 — kr^) ^dr^ + r^dQ 


(17) 


the Einstein equations (fT^ are reduced to the system 

. 2 j 

3 ^ — = SnGp + A, (18) 

P = -S% + P). (19) 

a 

Eq. (fT9l) results from the continuity condition Tl^.^ = 0, the dot denotes the time derivative, here 
and below the speed of light c = 1. 

Eor the ACDM with dust-like matter {p = 0) we use the solution of Eq. (fT^ p/pQ = (a/oo)”^ 
and rewrite Eq. (jlSp in the form 


a 


2 IJ2 
-^0 




( 20 ) 


Here the present time fractions of matter, dark energy (A term) and curvature 


^ 87rG/9(to) 

are connected by the eqnality 


Ha 


A 


H/j — 


k 

alHl 


T T f; 


resulting from Eq. ()20p if we fix t = to- 

If we introduce dimensionless time r and logarithm of the scale factor 


T = Hot, A = log ■ 


ao 


( 21 ) 


( 22 ) 


(23) 


equation (I^Up will take the form ^ -|- Ha + more convenient for nnmerical 

solving with the initial condition at the present time A\^^i = 0 equivalent to a(to) = ao- Here and 
below the present time t = to corresponds to r = 1. 

If we fix all model parameters, we can solve numerically the Cauchy problem for Eq. (1201) and 
calculate the values a(t)/ao, H{z), Dl{z) ([I]), dz{z) and A{z) ([7]). To compare them with the 
observational data from Ref. and tables niiHii we use the functions ([3]), (1131) and Cl and 
(under assumption about their Gaussian nature) the correspondent summarized function 


As — aIat + A/f + Xbao- (24) 

When we apply the ACDM model for describing the observational data from Sect. HI] (for 
z < 2.36), we use three free model parameters Hq, H^ and Ha (or H^, instead of Ha) and the 
additional parameter 


Pb{to) SirCpbito) 

ilh - - - - 7^ - 

Per 3F7| 
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if we use the fitting formula (IlOp for rs{zd)- For the formula (|lip the baryonic fraction is not a 
model parameter, we mentioned above that the value r^-Zi plays the role of an additional parameter 
in the case (jlll) . In both approaches we test the ACDM model and estimate the best ACDM fit 
correspondingly for Qb and • h. 

The results are presented in hgure [H where dependence of the best (minimal) value of the 
function (I24p minx^ = min xi; is shown in the top-left panel for the case (llOh . Here 

[—I |—I 

and for all models we assume '^rriy = 0.06 eV [71, Isj. One should note that the formula (fTop is 


m- 


insensitive to a neutrino contribution in the range — 0-23 eV 

For the variant with Eq. (jlip the similar dependence of min xf; on • h is in the top-right panel. 
Here and below we draw these graphs as solid red lines for the fitting formula (|10p and as dashed 
blue lines for the variant (HU). 



FIG. 1: The ACDM model: dependence of the best fit min Xs on Hb for the fitting formula (fTOll (the 

Hq ,r2m 

left panels) and on rd ■ h for Eq. (Ilip (the right panels). In the bottom panels we show the correspondent 
dependencies for parameters Hq, 17^, Ha, Hfc of the Xs minimum point. 


The bottom panels of figure [T] illustrate how coordinates h, Qm, Ha and = 1 — ~ Ha of 

the minimum point for min Xv depend on the correspondent parameters H^ and ■ h. 

One can see that for the formula (jlOp dependence of minxs on Qb is rather sharp: we have the 
distinct minimum at the value (1121) Qb = 0.044. Below we use this value as the fiducial one for all 
models. The correspondent dependence in the top-right panel is more smooth, however it results 
in the optimal (fiducial) value r^ - h = 104.57 ± 1.44 Mpc in the formula (jlll) . 

If we fix these parameters as described above, we can test the ACDM model for different values 
of the remaining 3 parameters: Hq, and Ha. The results of calculations are presented in 
tables lYllYllYl] and in figure [2j In the top-left panel of the figure we see how minimum of the 
function (|24l) minxs = As(-f^o) depend on the Hubble constant Hq: red solid lines in the 

top panels describe the model with the formula (HOI) , blue dashed lines correspond to the variant 
(HU- For the sake of comparison we present here graphs for the fixed value rs(zd) ([9]) as green lines 
with dots. These minima are calculated for each fixed value Hq. 

One may see that the function minxs(77o) for the fitting formula (jlOl) (the red line) has the 
maximal spread and achieves its minimum minxs — 577.39 at 

-1 


Hq = 70.07 ± 1.82 km (s Mpc) 


(26) 
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FIG. 2: For the ACDM in the top panels we present how minx|; depends on Hq, flm, and 0^ with 
different fitting formulae for for Eq. (fTOll with fit, = 0.044 (red solid lines), for Eq. ([TT1) (blue dashed 
lines), for Eq. (|9|) (green lines with dots). The correspondent dependencies for parameters of Xs minimum 
point for Eq. m are shown in the 2-nd row. In other panels Icr, 2a and 3a level lines are drawn for xh 
as red filled contours for Eq. (fTUl) and blue filled contours for Eq. m and also as black lines for Xsn 
green lines for x'sn + X^i optimal values of the third parameter are fixed and shown. 


for the following values of other parameters: Qm — 0.278, ~ 0.75 (presented in table HVl) . It is 

interesting to see in table HVl and in figure [2] that for the fitting formula dill) (blue dashed lines) the 
minimum min^s the optimal values of all parameters are the same, there is some difference 
only in Icr errors. To estimate Icr errors in the equality d26p we use the one-dimensional likelihood 
function Cy,{Hq) ~ exp(—x|/2) corresponding to min x|(i7o). 

Optimal values and errors for other model parameters are calculated similarly, they 

are tabulated below in table ITVl in comparison with other estimates and also in the next section in 

























































































































10 


tables |V] and |Vl] correspondingly for the expressions (fTOll and (fTTI) . 

Estimations of and in tables [TYIEIEI] are connected with graphs in the next two 
panels in the top line of figure [2l which present how minima of Xe depend on and on fij\. 
In particular, the solid red line in the second top panel describes miny; 2 (nm) = niin Xe'j the 

corresponding likelihood function CE{^m) ~ exp(—x|;/2) determines the la error — 0.008 in 
tables imiYi This panel demonstrates that dependencies of min Xe on for different variants of 
rs(zii) are rather close and have rather sharp form. It is connected with the contribution in Xbao 
from the value A{z) ([7|), because A{z) is proportional to and Xbao sensitive to Qm 

values. For the red line in this panel we have the additional dependence on in the formula 
(jlOp . In the third panel the functions minxs(f^A) with distinct minima have some difference in 
their the la errors. 

Estimations of in tables m-M are calculated via the function minxs(f^fc) = niin Xe- 

Ho,Qrn 

These graphs are shown in the top-right panel of figure [2j 

The panels in the second row in figure [2] demonstrate how parameters of a minimum point of 
the function Xe with the fitting formula (fTO]) depend on Hq, Hm, ^k- In the left panel these 
coordinates are Qm and Haj but also the value = 1 — — Ha is drawn as the black dashed line. 

The optimal value remains practically constant in contrast with Ha and H^. The graphs of 
/i(Hm),..., h(Hfc) in other panels show h = Hq/100, where Hq is the optimal value corresponding 
to the minimum point of Xe- 

In other panels of figure[2]we present the results of calculations as level lines at la (68.27%), 2a 
(95.45%) and 3a (99.73%) confidence levels for the functions x‘^{Pi^P 2 ) in planes of two parameters, 
if the third parameter is fixed. For example, the functions x|;(Hm, Ha) for the fixed optimal value 
(j26p of Hq are shown in the bottom-left panel of figure [2] as red hlled contours for the formula (jlOh 
and blue filled contours for Eq. (HU). The corresponding level lines for x‘sn for Xsn + Xh 
shown as black and green lines. 

In other panels only Xe liUsd contours for the cases (lIOp and (llip are compared for different 
pairs of parameters. Points of minima for the functions Xe marked as red (or blue) circles, we 
mentioned above, that they coincide for the variants (jlOp and ()II|1 . In two bottom-left panels the 
difference is in the fixed parameter (Ha or H^.) and in a choice of the foreground between the cases 
(fTOl) and (fTT]l . 

Our estimations of the ACDM parameters for two variants (jlOp and (llip of the fitting formula 
for Trf are to be compared with the the following best fits for these model parameters from surveys 
of the Wilkinson Microwave Anisotropy Probe (WMAP) and Planck Collaboration 0,0 in 
table HVl 



This paper 

WMAP, Planck surveys 


Td (Ho]) 

Td (HU 

WMAP 9y [8] 

Planck 13 [9] 

Planck 15 [Ihl 

Ho 

70.07 ± 1.82 

70.07 ± 1.27 

69.7 ±2.4 

67.3 ± 1.2 

67.8 ±0.9 


0.278 ±0.008 

0.278 ±0.009 

0.279 ±0.025 

0.314 ±0.02 

0.308 ±0.012 

Ha 

0-750t°:°^i 

0.750 ±0.034 

0.721 ±0.025 

0.686 ±0.025 

0.692 ±0.012 


-0.028t°:!]l° 

-0.028t°:!]i 

—0 nr)O7+0-0039 

U.UUZ/Q QQ3g 

-0.0005t°;™65 

-0.005t°:°)? 

H{, 

0.044 ±0.004 

- 

U.U^OC)_0.0024 

0.0487tHm? 

0.0484tO;°°)4 


TABLE IV: Estimations of the ACDM parameters. 


One can also add the estimates for the fixed ([9]): Hq = 69.27 ± 0.93 kms“^Mpc“^, H^ = 
0.280 ± 0.009, Ha = 0.732lgQ55, H^ = —0.012 ± 0.045 (corresponding to green dots in the top 
panels in figure [2]). In the case ([9|) the la error for Ho is smallest, because depends on h. 
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We see that our estimations of the model parameters in table IIVI for the cases (1101) and (|lip are 
in good agreement with the WMAP estimates, but they are in la or 2a tension with the values 
of Planck Collaboration. On the other hand, all ffo values in table m have essential tension with 
the Hubble Space Telescope group 581] estimation: Hn = 73.8 ± 2.4 kms“^Mpc“^. 


The latter value was used as a prior in Refs. [55|, [5^ for describing Type la SNe, BAO and H{z) 
data with the help of the ACDM, XCDM and (/)CDM models. One may conclude, that for the 
ACDM this choice of Hq was unsuccessful in comparison with another value Hq = 68 kms“^Mpc“^, 


chosen in Refs. 55|, |56| 


IV. MODIFIED CHAPLYGIN GAS 


In the model with modified Chaplygin gas (MCG) this gas has the following equation of state 

iil-Iij 


Pg = wopg - B p~ 


(27) 


for its density pg as a part in the total density (I16|) . MCG can unify dark matter and dark energy. 
If rco = 0 the MCG model with Eq. (12711 is reduced to the model with generalized Ghaplygin gas 
(GCG) with EoS 0, 


Pg = -Bp-^. (28) 

In our paper the GGG model was applied to describing the observational data for Type la 
supernovae [^, H{z) with 34 data points and BAO with 7 data points for dz{z). In this paper we 
consider the enlarged number of data points from tables mi mil the more general MCG model (1271) 
(in comparison with the GCG case wq = 0) and also we calculate the function Xsn ® with the 
covariance matrix Csn- 

The MCG and GCG models are to be explored as two-component models with usual dust-like 
baryonic matter component pb and the Chaplygin gas component pg with EoS (j27)) . In this case 
the total density (fT6]) is 


p = Pb + Pg, Pb = 0. 


(29) 


However the first component pb and the corresponding fraction ()25p may include not only visible 
baryonic matter but also a part of cold dark matter with p = pdm- Our practical applications of 
these models in Ref. [1^ and in this paper (see the top-right panel of figure [3|) demonstrate rather 
weak dependence of minxs on Hf, for the model assumption (ITT]) , but the strong fife dependence for 
the fitting formula (jlOp . This behavior resembles the ACDM model, where separation of baryonic 
and cold dark matter in appears only in Eq. (llOp . In the MCG and GCG models these matter 
fractions may also be mixed in their observational manifestations, so below we did not use Dfe as 
an usual free parameter of the theory, but fix its fiducial value (I12p in the main part this research 
(except for calculations, presented in right panels if figures [3] and d]) . 

Equation (fl^ for the MCG model (fTTl) is integrable, so the analog of Eq. (IMl) for this model is 




^ = nb(-) ' + (i-Dfe-Dfc) 

Here the dimensionless parameter Bg = B Pq^~°‘ / {1 -|- rco) is used instead of B. Thus in the MGG 
model ((271) we have 6 independent parameters: Hq, Dfe, Qk, wo, Oi and Bg- Naturally, the GGG 
model has 5 parameters: the same set without wq. 


^ao 




(30) 
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In the MCG model the parameter 0.^ from its formal definition (|21l) equals = 1 — flfc in 
accordance with Ha = 0 and Eq. (j22p . However, the expression A(z) ([7]) has the factor \/n^, so 
one should use the effective value in any model. If we compare the early universe limit z ^ 1 
in the MCG equation ([30]l with the ACDM equation (f20]l . we obtain the effective value 

+ (1 - (31) 

But for the majority of observational data in Ref. and tables miiiiii redshifts are 0 < z < 1, 
so to describe correctly these data we have to consider the present time limit of Eq. (j30|) . If we 
compare limits of the right hand sides of Eqs. (I20|) and (1301) at z —)• 0 or M 0, we obtain another 
effective value 

+ (1 - Hb - Hfc)(l - Bs){l + wo). (32) 

Expressions (|3T]l and (f3^ and their contributions in Xbao As were compared in Ref. . 
In this paper we use below Eq. (|32]) and its analogs for other models. 

Figure [3] shows how the MCG model (|27n describes the observational data from Ref. and 
tables im mil in comparison with the GCG model (j28p . Notations are the same as in figure El In 
the top line panels red solid lines and blue dashed lines denote graphs of miny;|. for the MCG 
model with expressions (llOp and (jllli correspondingly. The similar functions for the GGG model 
are shown as orange dash-dotted lines, if rs{zd) is calculated with the formula (fTOj) and as violet 
dashed lines for the case m- 

These graphs in the top line describe how Xs mentioned models depends on one chosen 

model parameter: Hq, wq, a, Hf,. All these curves determine the optimal values and errors of 
the model parameters in tables |Vl IVII The second row panels of figure [3] correspond to panels in 
the top line and present dependencies of coordinates of minima points on Hq, ..., for Xs fo 
MGG model with the fitting formula (fTO]) . 

In the top-right panel of figure [3] the value min x|; for the MCG model means the minimum 
over other 5 parameters for each fixed H;,: minxs = ™fo As- Our calculations support 

_ Ho,Q.i^,WQ,a,Bs 

the previous conclusion [39| about weak dependence of minxs on for the case (fTTI) with = 
{f'dh)fid ■ h~^ as for the the fixed value ([9]) for both models. This is connected with possible 
mixing of baryonic and cold dark matter in Hft. However, for and the fitting formula (lIOp we see 
the sharp dependence of minxs on Qb] for both MGG and GGG models this picture coincides with 
behavior of the AGDM model in figure [H As mentioned above, we will not use Qb as an usual free 
model parameter and we fix its ACDM fiducial value (1121) . It is interesting that for the MCG and 
GGG models this value Qb — 0.044 ± 0.004 is very close and practically coincides with the AGDM 
fiducial value of this parameter. 

In figure [3] the parameter Qb is varied only in 3 right panels, in other 12 panels this value is 
fixed in the form (II2p . In particular, the value minxs(77o) for the MCG model in the top-left 
panel is the minimum over 4 parameters min Xs with fixed Qb = 0.044. The similar picture 

Uf;,wo,a,Bs 

takes place, when we study dependence of these minimal functions on wq, a, for example, 
minx|](Hfc) = min Xs- For the GCG model with wq = 0 these minima (orange dash-dotted 

Ho,wo,a,Bs 

lines and violet dashed lines) are taken over 3 remaining parameters. 

The graphs of minxs(77o) fo the top-left panel with two expressions for rs{zd) resemble the 
AGDM case in figure El for both MCG and GCG models. For all three considered models this curve 
for the case (jllj) is slightly more sharp than for the formula (lIOp . this behavior may be seen in 
tables lYllYll where all models are compared. Unlike the ACDM model, for the MCG and GCG 
models we have different optimal values of Hq (and also of a and other parameters) in the 
cases m and m- 
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FIG. 3: For the MCG model dependencies of minx| on one chosen parameter iJo, Hfc, wo, ct, fib are 
presented for the fitting formulae (ITU)) (red solid lines) and (fTT)) (blue dashed lines) in comparison with the 
GCG model (orange dash-dotted lines and violet dashed lines for these two choices of r^). In the second 
row there are correspondent MGG parameters of a minimum point for with Eq. (1101) . In other panels for 
the MCG model Icr, 2cr, Scr level lines are drawn in the planes of 2 parameters (other parameters are fixed). 
Notations correspond to figure [2j 


Obviously that for the GCG model the graphs of Xe in the top panels of figure [3] lie higher, 
then the correspondent MCG graphs. These curves converge at points, where the optimal value of 
tco equals zero in the panel below. 

The graphs of minxs(IlA:) in the second top panel are asymmetric, for > 0 these values for 
the GCG and MCG models are close, but for 0^ < 0 the curves for these models diverge for both 
variants of rs(zd)- In the third top panel we have asymmetric dependence of these minima on wq 
for the MCG model. 

Dependence od minx|j on a is essentially different for the considered models: for the GCG 
model these curves have sharp minima, in particular, for the case (uni) at a ~ —O.I with a ~ O.I, 
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but for the MCG model the dependence is rather smooth, the minimum at a ~ 0.54 has a ~ 0.9. 
These results are also presented in table correspondent values for the function m are in 
table ED One may conclude that changes of the parameter wq may compensate changes of a. 

In 5 bottom panels of figure [3] for the MCG model we present la, 2a and 3a level lines in the 
planes of 2 model parameters (in notations of figure [ 2 |), in particular, red filled contours denote 
levels of with the fitting formula (llOj) . For this case in the bottom-left panel in the {a, Bg) plane 
contours for XsN (black) and Xsn + xjj (green lines) are also shown. Contours for with the 
formula (jlll) in all panels are drawn as blue lines. Other model parameters are fixed and shown in 
the panels and in table El they are optimal for the fitting formula (jinj) (but not optimal for the 
case (fTTTl l. 


Model 

minxi 

Ho 


other parameters 

ACDM 

577.39 

70.07 ± 1.82 

-0.028^^48 

Qm= 0.278 ± 0.008, Ha = 0.750lo;o55 

GCG 

576.13 

69.46 ± 1.88 

0.0261°:°^^ 

a = -0.100/:°:°^°, Bg = 0.738t°;“4 

MCG 

575.60 

70.19 ±2.15 

-0.044 ±0.122 

a = 0.538;°J°2, Bg = 0.714/° 
wo = -0.166/°;322 

EoS 

(IHl) 

575.15 

70.28/|?2 

— 0 04 

U.U^O_Q Qg4 

Po = -0.904/0;256, p = -0.042/°:°4i, 

wo =0.183 ±0.222 


TABLE V: Models and la estimates of model parameters, if rs{zd) has the form (flOl) : Hf, = 0.044. 


Model 

min x| 

Ho 


other parameters 

ACDM 

577.39 

70.07 ± 1.27 

-o.o28/!]:011 

^m= 0.278 ± 0.009, Da = 0.750/[];J]|^ 

GCG 

576.48 

70.27 ± 1.28 

-0.009 ±0.040 

(A = -0.069 ± 0.074, Bg = 0.753 ± 0.013 

MCG 

575.61 

70.44 ± 1.30 

-o.o6o/°:°?5 

a = 0.613/°;«°4, Bg = 0.716/°;°^, 
wo = -0.176/0//^ 

EoS 

(|M|1 

575.14 

70.42 ± 1.30 

-0.048/n48 

PO = -0.928;°|42^ p = -0.045/rof9, 
wo = 0.205/°;!®^ 


TABLE VI: Models and Icr estimates of model parameters with rs{zd) from Eq. (jUJ; fib = 0.044. 


For the GCG model (|28p the values in table ED are close to our previous estimations 


--1-0.072 


Bg 


= 0.759^Q;Qlg with 7 and 34 data 


Ho = 70.093 ±0.369, Qk = -0.019 ± 0.045, a = -O.O 66 Io. 074 , 
points for dz{z) and H{z) correspondingly. 

These estimates for the MCG model should be compared with similar results for this model in 
The authors of Ref. 451 for the flat model with = 0 described the observational 


papers 


data with 557, 15 and 2 data points correspondingly for supernovae, H{z) and BAG, but they also 
included the cluster X-ray gas mass fraction data. Their la estimations Hq = 70.711 


+4.188 


narrow box ^ = 0.00189lQQQ75g lie inside our la estimates. 


-3.142 

Bg = 0 . 77881 q;q 723 more wide then ours in tables El ED however a = 0.1079^q; 2539 and the 


and 


In Refs. [46], l47|] the MGG model with flfc = 0 is applied for describing 12 H{z) data points, 11 


points for the growth function / = d log (5/d log a of the large scale structures, 17 points for as{z) 
and also in Ref. 471] observations of supernovae and 1 data point for BAG. The authors did not 
demonstrate their estimates for Hq, but noted ambiguously function for the background test 
is minimized by the present Hubble value predicted by WMAP7”. The best fit values of other 


parameters in Ref. [47|] are wo = 0.005, a = 0.19, Bg = 0.825 with errors, calculated for pairs of 


these parameters. Gnly the estimate for Bg is in tension with our results in tables El ED 
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V. MODEL WITH QUADRATIC EQUATION OF STATE 


It is interesting to compare the MCG model and the model with quadratic equation of state 




Pg=PQ+ WoPg + j3pl, 


(33) 


because both models have 6 parameters: Hq, Qfc, and 3 parameters in the EoS (|27)) or (f33ll . It 
is convenient to use the critical density pcr = 3Hq/{8ttG), introduce the dimensionless parameters 
Po = Po/Per, /3 = (3Per instead oi po, f3 and rewrite Eq. (|33]l in the form 

Pg = POPer + WoPg + ^p\l Per- (34) 

Similarly to the GCG and MCG models the model with quadratic EoS ()33|) or (IM)) has the 
analytical general solution of Eq. (I19p 


P^ 

Per 


= < 


^ r r-^tan(§^.4) _ , _ , ' 
2/3 _l+|A|-i/2tan(|vlA|.4) ° 

^ (|^+ l/r)“^ - I -Wo 

p-{Q,m-p+){a/ao)~^'^-p+{^m-p-) 

(Om—p+)(a/ao)“®'^ —f 2 m+P- 


, A<0, 
A = 0, 
A > 0, 


(35) 


depending on a sign of the discriminant A = (1 + rco)^ — 4/3po- Here 


— 1 E — ‘213^'fYi “b 1 “b Wo, p± — 

The equation (fT8]l for this model is reduced to the form 

dA-' 


— 1 — rco ± \/A 


wq 

2/3 


^ = log ■ 


ao 


dr J Per 


Hi 


(36) 


We solve this equation numerically from the present time initial condition A\^^i = 0 “to the past”. 
We can use analytical solution ([35]) or solve Eq. (fTOl) numerically, these approaches are equivalent. 
Eor the model ([3l]l the effective value 


-)- pq -)- Q,.^(I + uiQ + P^m) 


(37) 


is calculated in the z —>■ 0 or ^ 0 limit similarly to the MCG model. 

The model with quadratic EoS (1341) in the domain /3 > 0 may have the following singularity 
in the past: when t —)• + 0, density grows to infinity, but the scale factor remains finite and 
nonzero: lim p„ = 00 , lim a = a(tA A 0. This behavior resembles the Type III finite-time future 


singularity from the classification [361. l37[|. In Ref. [53(] the author did not see these singularities, 
because he considered only negative values /? = —(wo + l)/pp. 

In the bottom-right panel of figure S] we present the example of singular solution as the red 
dashed line for o(r)/ao and the black dash-dotted line for 0.01 • pir)!per- This singularity is 
compared with the regular solution (the blue solid line for a(r)/ao) with the optimal values of 
model parameters from table jVl Eor the singular solution in this panel /3 = 0.02, but other model 
parameters are from table jV] 

We have to exclude these nonphysical singular solutions, for this purpose we can use different 
approaches. The simplest way is to add the penalty contribution in in the form 


= -Pi[exp (P 2 a(t*)/ao) - !]■ 


(38) 
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The function (l38t) with Pi = 0.3 and P2 = 20 successfully helps to avoid this singularity. For 
regular solutions a(p) = 0 and the contribution (f38ll vanishes. In the top panels of figure |4] 
functions min Xs with the contribution (j38p are shown as solid red lines for the formula (jlOp and as 
blue dashed lines for = {rdh)fid ■ for the case without penalty (f38P these graphs are dotted 
lines of the correspondent color: red for the formula (IlOp and blue for (|lll) . 

More natural method to eliminate nonphysical singular solutions is to include early time param¬ 
eters (for example, Zd) into consideration. Unfortunately, both equations (fTOl) and ifTTI) or rs{zd) 
are insensitive to the mentioned singularities. So we can take into account the cosmic microwave 
background (CMB) constraints, in particular, for the values Q 

V = (wfe, Web, PM(1090)/rd), (39) 

where ui = Dm{z) = Dl{z)/{1 -|- z). If we calculate the vector (p^ . compare it with the 

estimation [3,1^ Av = v — (0.02259, 0.1354, 94.51) and add the corresponding term to the Xs 
function (I24p . we obtain 

x|+= x| + Av-C^|^5 (Av)^. (40) 

Here the covariance matrix from Ref. 0 contains cn = 2.864 • 10 C12 = —4.809 • 10 

ci3 = -I.Ill • I0-^ C22 = 1.908 • I0-^ C23 = -7.495 • 10"®, C33 = 0.02542. 

Graphs of one parameter functions minx|]_|_ with the CMB contribution for the case (llOh are 
also presented in the top panels of figure 0] as black dash-dotted lines. 

We demonstrate in figure 0] how the model (j34p is effective in describing the observational data 

from Ref. 0 and tables miiiiii Notations are the same as in figure [3l Here we also fix the value 

(fT^ = 0.044 (except for 2 top-right panels) and do not use Hfe as a fitting parameter. The 

dependence minxs(Hb) = min Xe (shown in the top-right panel) is rather weak for the 

/fo,Ofe,ioo,po,/3 

case (ffH) I’d = {fdh)fid ■ h but it is essential for the formula (flOP and for the function minx|]_|_. 
In the case minxl with Eq. (|I0I1 (the red curve) the dependence minxKUb) is very close to the 
correspondents functions for the ACDM, GCG, MCG models in figures 01 [3l its minimum is also 
at Qb - 0.044 ± 0.004. 

In the top line panels of figure 0] we draw graphs of Xs minima depending on one model pa¬ 
rameter {Hq, Hfc ...) for the model (l34)l with the penalty function ([38l) as solid red lines for the 
formula m and as blue dashed lines for the case dm)- Correspondent lines without contribu¬ 
tion ()38p are shown as dots of the same color. In other words, if we calculate minxs only for 
(physical) regular solutions, we obtain the solid or dashed lines; for dots we also include singular 
solutions without physical interpretation. These lines coincide in domains where best values of /3 
are negative and corresponding solutions are regular. One can see in figure 0] that optimal values 
of model parameters correspond to regular solutions in both considered cases (fTOl) (table IVj) and 
fn]) (table EH). 

Minima of Xe Xs+ here have the same sense as in figure[3]for the MCG model, in particular, 
in the top-left panel minxE(F7o) = min Xs- Dependence of this minimum on Hq for two 

formulas for rs(zd) resembles other considered models (compare with figures [21 [3]) : the graph for 
the case (|III) is more sharp, than for Eq. cnD, it corresponds to larger optimal value of Hq and 
smaller la error in table IVII than the values in table El For other panels of the top line in figure 0] 
dependence of minxs on parameters Hfc, wq and po results in the correspondent the Icr estimates 
in tables El ED 

The functions min Xe+ with the CMB contribution (1401) depending on one parameter are shown 
as black dash-dotted lines in the top line panels. We see that the absolute minimum of this function 
is minxE+ — 576.16, in a bit exceeds the corresponding value minxs — 575.15 for the case (|I0l) 


17 



FIG. 4: For the model with quadratic EoS (j34|) one parameter dependencies of minxl with rs{zd) in the 
forms (flUl) (red solid lines) and (fTTl) (blue dashed lines), of minx|, with CMB contribution and Eq. (flUl) 
(black dash-dotted lines) and also coordinates of minima points, level lines are presented in notations of 
figure [3l In the bottom-right panel scale factors a/o(0) for the regular solution (the blue solid line) and the 
singular solution (the red dashed line) are shown. 


in table |Vl However, the CMB contribution in the form ()40p works as a very narrow filter for 
some model parameters, in particular, for the considered model with EoS (13411 this contribution 
rigidly constrains the value flfc, so we have the Xe+ estimation = —0.090 ± O.OOI. The curve, 
corresponding to this narrow range, is shown in the second top panel of figure 01 it is too narrow, 
so it looks like a vertical black segment. 

The x|_|_ estimations of other parameters are also more narrow (see the top line in figure0]), than 
for both cases of Xe’ connected with restrictions for Q/. and other correspondent parameters 
for other models. Some of these restrictions look like artificial od connected with the concrete 
choice of the CMB vector (|39p . So further we consider the functions miny;|; without the CMB 
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contribution (j40p with the fitting formula (1101) as the most reliable indicator. 

Panels in the second row of figure |4] correspond the upper panels for the function min Xs with the 
fitting formula (jlOp with the penalty contribution (|38l) . They demonstrate evolution of coordinates 
of the minimal point if we vary the chosen parameter. In 4 bottom panels functions depend on 
two model parameters (p0)/3; Po-,Hq\ when other parameters are fixed with their 

optimal values from tabled Red filled contours denote the case (fToP . blue lines corresponds to 
(HH). These level lines demonstrate similarity with the MCG model in figure El but for the model 
(I34p we have alternative behavior in the plane. One should emphasize that all lu, 2a and 

3a level lines for in figure 0] and optimal values of the model parameters in table |Vl IVII lie in 
the domain with regular solutions of the quadratic EoS model (p^ . 


VI. CONCLUSION 


In this paper the Type la supernovae observational data from Ref. and estimations of BAO 
parameters and H{z) from tables Ull HIl are described with the models ACDM, GCG ([28]), MCG 
(127h and the model with quadratic EoS (I34j) . Two approaches in estimation of the sound horizon 
scale rs{zd) are used and compared: the fitting formula ([TUP with results, tabulated in table IVl 
and the expression (fTTp = {rdh)fid ■ h~^ (table IVTl Optimal values of model parameters with la 
errors in these tables are calculated via one-parameter distributions (figures [2]-|3]) . 

We also considered the CMB contribution in the form (|40h , the results for the model with EoS 
()34l) are shown in figure 01 However, this contribution appeared to be too sensitive and restrictive 
for some model parameters, in particular, (or Ha for the ACDM model). So we consider the 
estimations of model parameters in tables |V] and IVII as more reliable. 

It is interesting that predictions of the ACDM, GCG, MCG models and the model with quadratic 
EoS for flf, are very close (fib = 0.044 ± 0.004), if we adopt the fitting formula (fTOl) for rs(zd)- We 
use this fact and do not consider fib as an usual model parameter and fix its value in the form 
m- One should note, that predictions of different models for Hq and flk in tables |V] and IVj are 
also rather close. 

Absolute minima of Xe with the formula (llOj) in table |V] differ from the correspondent minima 
of in table rvTI with Eq. (Hip , but the hierarchy of all considered models is the same in these tables. 
In particular, the absolute minimum of x^ in table 13 vary from the worst value 577.39 for the 
AGDM to the best result 575.15 for the model with quadratic EoS ([34P . Note that the advantage 
of the MGG model in comparison with GCG is larger in the case (fTTI) in table IVTl 

However, effectiveness of a model essentially depends on its number Np of model parameters 
(degrees of freedom). This number is used in model selection statistics, in particular, in the 
following Akaike information criterion 57|, 


AIC = minxi; + 2Np. 

If we fix the value fib in the form (|12l) for the models GGG, MCG and with EoS (|34p and do 
not use this parameter as a degree of freedom, we will have the numbers Np and AIC for Xs from 
table 13 for the considered models tabulated here in table IVHl 

This information criterion works against models with large Np and adds arguments in favor of 
the ACDM model. 
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Model 

min 

Np 

AIC 

ACDM 

577.39 

3 

583.39 

GCG 

576.13 

4 

584.13 

MCG 

575.60 

5 

585.60 

EoS dSH) 

575.15 

5 

585.15 


TABLE VII: Akaike information criterion for the models. 
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